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Abstract 

We extend recent low temperature analyses of competing orders in the cuprate superconductors 
to the pseudogap regime where all orders are fluctuating. A universal continuum limit of a classical 
Ginzburg-Landau functional is used to characterize fluctuations of the superconducting order: this 
describes the crossover from Gaussian fluctuations at high temperatures to the vortex-binding 
physics near the onset of global phase coherence. These fluctuations induce affiliated corrections 
in the correlations of other orders, and in particular, in the different realizations of charge order. 
Implications for scanning tunnelling spectroscopy and neutron scattering experiments are noted: 
there may be a regime of temperatures near the onset of superconductivity where the charge order 
is enhanced with increasing temperatures. 
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I. INTRODUCTION 



A number of recent perspectives El El H , 



l a LZ | have h ighl ighted new exp erimental 
B El El El El El and theoretical )Mlm O LlJllI, H H B 0, HI M, works 
exploring the interplay between the multiple order parameters which characterize the ground 
state of some of the cuprate superconductors. Good evidence was obtained for a strong 
coupling between the superconducting order and density wave order in spin/charge/bond 
correlations (described more precisely below). In particular, by tuning the superconducting 
order by an applied magnetic field at very low temperatures (T), a strong field- dependent 
variation was observed in the latter correlations. 

In this paper, we explore the possibility of observing related connections in the finite 
temperature 'pseudogap' region above the superconducting critical temperature, T c . Here, 
the superconducting order has strong T dependent fluctuations; we will compute these 
fluctuations in the framework of a two-dimensional Ginzburg-Landau theory, including a 
precise characterization of strong fluctuations obtained from numerical studies. We will 
show that the model of Ref. [ly predicts that such fluctuations lead to a corresponding 
sympathetic variation in the autocorrelations of the other orders. Working to linear order 
in the coupling between superconductivity and these orders, we provide a computation of 
certain universal characteristics of the T dependence of the latter fluctuations. Our results 
will also be formally extended to T < T c for completeness, but it must be noted that we 
neglect the inter-layer coupling and quantum effects, which become important at lower T 

We begin by defining the order parameters under consideration. The primary order 
is the complex superconducting order \I/(r) which describes the spatial variation in the 
order associated with condensation of Cooper pairs. This is expected to undergo strong 
'phase' fluctuations [29( for T near T c . Using t he p roximity of the underdoped cuprates to 



31| argued that 'amplitude' fluctuations 



a superffuid-insulator quantum transition, Refs. |3fJ, 

should be treated at an equal footing (32| , and proposed that such thermal fluctuations could 
be described by a classical partition function of a suitable universal continuum limit of the 
Ginzburg-Landau free energy: this will be reviewed here in Section ITTT1 Such an approach 
describes the crossover from Gaussian superconducting fluctuations at temperatures well 
above T c , to the vortex physics of the Kosterlitz-Thouless transition near T c . A dynamic 
theory with a similar static component (although with a lattice cutoff) was recently used 
jH Q,!!^ to describe the notable measurements [36| of the Nernst effect. 

This fluctuating superconductor is also expected to have appreciable correlations in other 
order parameters. The spin-density-wave order is described by the complex 3-component 
vectors <3> yQ ,, where a = x, y, z extends over the 3 spin directions, and the spin operator 
on site r, S a (r) is given by 



S n 



Re 



r <Mr) + e l ^ r <M r ) 



Here K SXj3/ are the spin-density-wave ordering wavevectors along the x and y principle axes 



(37r/4,7r) and K 



of the square lattice: near a doping of 5 = 1/8, we have K. sx 
(7r, 37r/4). In a similar manner we can define bond order parameters 
the modulations in the exchange energy of a pair of spins separated by a distance a: 



*,</ 



S a (r)S a (r + a) = Re [e lK -' r ax (r) + e iK ^ r 



by examining 



(1.2) 
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The special case a = of 4>ea,y is a measure of the charge density wave order. Comparison 
between ()1.2|) and suggests that the ordering wavevectors are related by K. CXjV = 2K. sx>y , 
and this is observed experimentally. 

A number of other order parameters which are invariant under spin rotations, like <frax,yi 
can also be defined l2l| . These include the site charge density, the average electron 



kinetic energy in a bond, or modulations in the pairing amplitude. By symmetry, all such 
quantities will have modulations at the wavevectors K CXiy , and we can therefore expect 
that their order parameter fluctuations will track those of <fi ax , y . Differences in microscopic 
physics can, of course, make some of these modulations much larger than others. We will 
not explicitly consider all such possibilities here, and the reader should view 4> a x, y as a suit- 
able representative of the order parameters characterizing modulations at the wavevectors 
K. cx y in observables invariant under spin rotations. We will subsequently refer to the order 
represented by simply as charge order. 

While the focus of this paper is on the interplay between the superconducting fluctuations 
and the orders mentioned above, it should be clear to the reader that our considerations are 
quite general. Simple extensions lead to similar effects in the interplay of superconductivity 
with most of the other orders in the zoo of possibilities considered in the theory of the 
cuprates. 

In considering correlations of $ and in the fluctuation region, it is important to con- 
sider the influence of random static impurities which are invariably present in the cuprates. 
As almost all impurities preserve electron number and spin rotation invariance, their influ- 
ence on \1/ and $ will consist of perturbations in the random exchange class (this is discussed 
more explicitly in Section [nj. In contrast, the order breaks only lattice symmetries, and 
is consequently subject to the far more disruptive random field perturbations j3?|. In two 
spatial dimensions, this implies that true long-range order cannot develop as T — > 0, and 
that the correlation length saturates at a finite value. We will assume here that there is 
a local onset of 0, $, and \I> orders at temperatures where the pseudogap develops, but at 
lower temperatures correlations are predominantly controlled by the random-field disor- 
der, and have only a weak, intrinsic T dependence. This is also consonant with the result 
that thermal fluctuations are irrelevant at the random-field transition in higher dimensions 
(3?}. In contrast, the fluctuations of \I> and $ are strongly T dependent, and can have an 
infinite correlation length as T — > 0. The $ order becomes quasi- long-ranged at T = T c 
and has the strongly T-dependent Gaussian-to-vortex crossover noted above at T > T c . The 
$ order can also have the exponential rapid T dependence associated with the breaking of 
0(3) spin rotation symmetry as T — > 0. 

This paper will consider the regime above T c where 

(*(r)) = ; (<tW(r)> = ; (<Wr)> + (1.3) 

The non-zero value (0) is due to the presence of random-field perturbations which explicitly 
break lattice symmetries, and so allow to locally have a non-zero mean value which will 
fluctuate randomly as a function of r. As noted above, we assume that (0) only has a 
weak intrinsic T dependence. However, the fluctuations of the $, and orders are 
not independent, and so the strong T dependence associated with the Gaussian-to- vortex 
crossover in \1/ will induce a corresponding T-dependent variation in (0). This paper will 
compute this variation and suggest associated experimental tests. Strictly speaking, because 
there is only quasi-long-range order in \1/ below T c , the expectation values (jl.3|) apply also 
for T < T c : indeed, our methods and results extend also to T < T c . However, as noted 
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earlier, we neglect the effects of inter-layer couplings and of quantum fluctuations, and so 
our low T results should be treated with caution. 

Our theory for the fluctuating orders and their interplay is summarized in Section m 
which also contains are main results. Details of the continuum theory of the superconducting 
fluctuations and its Gaussian-to- vortex crossover appear in Section ITTT1 Section HV1 discusses 
experimental tests and possible extensions of our theory. 



II. CORRELATIONS BETWEEN FLUCTUATING ORDERS: MAIN RESULTS 

This section will introduce the free energies which control the fluctuations of the order 
parameters, and state our main results on the T dependence of the <fi order at T > T c . 

We describe the fluctuations of the superconducting order \l/(r) by a classical continuum 
partition function over the Ginzburg-Landau free energy [30| 

Z GL = J vm(v)e- rGL/{kBT) 

h 2 ,_ _ , , l0 , , , b , 



•F gl = I d z r 



V r ^(r)| 2 + a(T)|^(r)| 2 + -|^(r)| 4 



(2.1) 



We use here the notation of Refs |33|, |34j, |35f m*, a(T), b are parameters which can be 
computed, in principle, from the microscopic physics of the underlying electrons. The co- 
efficient of |\l/(r)| 2 , a(T), vanishes at a mean-field transition temperature, a(T c MF ) = 0, which 
will be distinct from the temperature T c at which there is a Kosterlitz-Thouless transition 
i.e. a(T c ) < 0. The purely two-dimensional, and classical theory (|2.1|) is expected to apply 
to the cuprates only for T > T c : below T c we have to also account for three-dimensional 
effects arising from inter-layer couplings, and for quantum effects at low enough T. All such 
effects will be neglected here, but for completeness, we will nevertheless discuss properties 
of the theory (|2.1j) over the full range of T values. 

An important point is that the functional integral in (|2.1|) is not defined on its own, and 
needs an ultraviolet regulator. In the physical system this is provided by the underlying 
electron physics on the lattice, but this is very difficult to characterize explicitly. Here, we 
shall follow the procedure proposed in Ref. [3(j: the ultraviolet dependence can be accounted 
for by a suitable renormalization in the value of a(T). However, because we do not know the 
explicit form of the ultraviolet cutoff, we cannot a priori compute the needed shift in a(T). 
This lack of knowledge can be circumvented by using the experimental value of T c as an 
input into our calculation. The knowledge of the actual T c , combined with the parameters 
in (J2.1j) then allows a quantitative computation of the Gaussian-to- vortex crossover with no 
free parameters. We re-iterate that (J2.1j) cannot be regarded as a fully predictive theory 
on its own, and so cannot, even in principle, predict the actual value of T c : once T c is 
determined by other means, precise quantitative predictions for other observables become 
possible. 

The Gaussian-to- vortex crossover can be expressed in terms of the following dimensionless 
parameter 



9 



h 2 \a(T) a{T c 



m*b 



k B T k B T c 

The parameter g should be a monotonically increasing function of T. For T T c , g ~ —1/T, 
at T = T c , we have g — 0, and above T c , g takes positive values. We will see later that 
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the present continuum theory eventually breaks down at large T, when g begins to acquire 
a non-monotonic dependence on T. The value of l/\g\ is a measure of the strength of 
corrections to the mean field theory of Zql- 

It is important to note that the T dependence of a(T) in (|2.1j) and (|2.2|) is non-universal, 
and this will lead to some non-universality in the T dependence of all our predictions. 
However, one of our main points is that there is a universal dependence on the parameter 
g. Moreover, once we assume a linear T dependence of a(T) near T c (as is commonly done, 
and we will do in (|3.13p ). the T dependence of our predictions becomes specific. 

Aided by the results of Ref. 13(1 l39l. we will show that it is possible to obtain precise 
predictions for a variety of correlators of Zql- We quote a result which will be useful in our 
analysis here of multiple order parameters: 



m* 



k B T c 



V(g,T/T c 



(2.3) 



where T>(g,T/T c ) is a universal function. The averages on the left-hand-side are evaluated 
under the partition function Zql at the indicated temperature. We will show in Section ITTTl 
that it is possible to re-express the two argument function T>(g,T/T c ) in terms of a single 
argument function F(Q) as in (|3.12|) . where Q depends upon g and T/T c as in (|3.11j) . Here 
we present results for the initial crossover from the Gaussian to the vortex regime, which 
occurs when g ^> 1: 



V(g,T/T c 



~ \n(380gT/T c ) + -i- \n(13.3gT/T c 
Z7T Z7r z g 



+ 



47r 3 g 2 



ln 2 (13.3( ? T/T c )-21n(7.86( ? T/T c ) + 0(1/ g 3 ). (2.4) 



The numerical constants appearing in the arguments of the logarithms are universal. These 
constants, and the constants appearing in the arguments of all subsequent logarithms, de- 
pend on only two universal numbers that have to be determined by computer simulations: 
the latter numbers are the constant Q c computed first in Ref. and the constant £ com- 
puted in Refs. [3^| H^. Additional higher order terms in ()2.4j) have also been computed 
and these will be presented in Section IIHI we show there that it is possible to account 
for all logarithmic terms that appear at higher orders in g. Numerical results for the full 
range of values of g appear in Section IIHI The expression ()2.3|) has ignored the possible 
T-dependencies of m* and b for simplicity: it is possible to account for these in a similar 
manner, as will become clear from the discussion in Section 11111 

It is worth noting here that vortices are already present in the Gaussian theory, associ- 
ated with zeros of \I/(r) jZJ]. The result (J2.4j) accounts for the initial correlations between 
these vortices, but does not include the vortex-binding physics of the Kosterlitz-Thouless 
transition. The latter is only accounted for by the numerical results in Section IIHI 

Let us turn now to the density wave order parameters $, <fi. The complete effective action 
for these order parameters has a rather complicated structure and was discussed in Ref. 
A simple Gaussian form will be satisfactory for our purposes here: 



d 2 r 



2 , t-2|$ 12 



+/Hte(r)^(p) + h% x {v)^{v) + (x - y) + 
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+ftyz( r )<M r ) + hl x {v)<f) ax {r) + (x ^ y) + . . .] (2.5) 



Apart from the usual Gaussian terms [16], the above contains complex random fields h$(r) 
and h$(r) which pin the 'sliding' mode of the charge density wave. These fields arise from 
impurities which preserve spin rotation invariance: as a consequence, notice that the random 
coupling is linear in the fields (f>, but that there is only a random-exchange coupling to 0(3) 
rotations in the spin density wave order. These simple facts have a number of interesting 
implications: 

(i) There can be no long range charge order in two spatial dimensions, even at T = 0. This 
implies that there can be no T = quantum critical point, tuned by the hole concentration, 
associated with the onset of such order. A quantum critical point associated with the 
restoration of 0(3) symmetry remains possible. 

(ii) The strong relevance of such random-field perturbations suggests that in the absence of 
couplings to other critical order parameters, the correlation length ^ can be assumed to be 
roughly temperature-independent at low temperatures. 

(in) The theories (|2.1j) and (|2.5j) describe a phase in which the expectation values in (jl.3|) 
hold. 

Finally, as promised, let us consider the influence of the ^ fluctuations described by Zql 
on the charge order correlations. The simplest coupling between the orders is a A | | 2 ( 1 aa: | 2 + 
|0aj 2 ) term, and, as in Ref. flH. this leads to the leading order correction 

&( T )=^o(T) + \(\*\ 2 ) T . (2.6) 

Here £</>o(T) is the 'bare' correlation length of the <fi order, which is expected to be only 
temperature dependent near T c . We input the value of (|^/| 2 ) as computed in ()2.3|) and 
Section IIII[ and obtain our main predictions for the superconducting fluctuation-induced 
modification in the correlation length. 



III. CONTINUUM THEORY OF THERMAL SUPERCONDUCTING FLUCTUA- 
TIONS 



This section will review the results of Ref |30| relevant to obtaining (J2.3)) and (|2.4j) and 
its extensions. Appendix [X] will review the work of Prokof 'ev, Ruebenacker, and Svistinov 
~39l 4oj ] on the dilute two-dimensional Bose gas and show that the results of their numerical 



simulations can be mapped onto universal quantities needed here. 

Ref. HO studied the following continuum theory of a N = 2 component real scalar ip a , 
a = 1,2: 



(3.1) 



(Here we have changed notation for the field, from $ a in Ref. |30|, to tp a here, to prevent 
confusion with the spin density wave order.) This theory maps onto (|2.1j) with the following 
corresp ondences : 



R 

U 



/ m*(ip 1 + itp 2 )/h 
2m*a(T)/h 2 
I2m* 2 b/h* 



(3.2) 
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It was argued 3(| that the continuum limit of required only the single renormalization 
of R to R: 

~ 2k B TU r A d 2 k 1 

R - R ~^^J a^WTr" (3 ' 3) 

Here we have introduced an ultraviolet cutoff A which is needed to regulate the theory 
T^. The renormalization in ()3.3|) is associated with logarithmic ultraviolet divergence of the 
one-loop 'tadpole' diagram; the renormalized R in the propagator on the right-hand-side 
accounts for tadpoles-on-tadpoles etc. All other diagrams are ultraviolet convergent, and 
hence the simple structure of the renormalization theory. 

It is important to note that (J3.3|) is the exact definition of R, and consequently R is not 
the fully 'self-energy' of the (p a field at zero external momentum; R only accounts for the 
resummation of tadpole graphs. In practice, the relationship (J3.3j) implies that, when we 
perform a Feynman graph expansion of any observable, we can ignore all tadpole graphs, 
and replace R by R in all propagators. Notice also that as the bare coupling R extends from 
— oo to oo, the renormalized coupling R extends from to oo. 

After the renormalization of R to R, all subsequent correlators of are ultra-violet 
convergent, and so we can safely take A — > oo in them. This implies that all these correlators 
are universal functions of the single dimensionless quantity that can be obtained from the 
parameters in (J3.1)) : this is the analog of the 'Ginzburg ratio', defined here as 

„ k B TU . . 

G = (3.4) 

For T <C T c , where i? < 0, we have Q — > oo. Conversely, for T > T c , i? > and £?—►(). 

The field theory (|3.1|) exhibits a Kosterlitz-Thouless transition at some critical tempera- 
ture, and the arguments above imply that this transition occurs at a universal critical value 
Q = Q c . The numerical studies of Ref. found Q c ~ 102. The value of Q c can also be 
obtained from the subsequent, and more precise, numerical simulations of Refs. 3^, 4(1 this 
connection is discussed in Appendix |XJ and (|A5J) yields 

£ c = 96.9±3 (3.5) 

We are interested here in the value of ((fi 2 )- This quantity requires a single additive 
renormalization before the continuum limit is obtained; hence we can write 

(ifl) f A d 2 k 1 . , 

" "' ° 1 +F{G) (3.6) 



k B T J 4vr 2 k 2 + R 

where F(Q) is a universal function. A number of analytic results for this universal function 
can be obtained from the methods of Refill, and details appear in Appendix |B| For Q — » 
(corresponding to T ^> T c ), perturbation theory in powers of U about the (p a = saddle 
point yields 

F(g^ 0) = (2.355711 x 1(T 4 )£ 2 + £>(£ 3 ) (3.7) 

All subsequent terms in the above expansion involve only integer powers of Q and there are 
no logarithms. For Q — ► oo (corresponding to T T c ), we expand about a saddle point 
with ip a 7^ 0. As shown in Ref. Hi} 1 , this is done by introducing a 'dual' coupling Q D related 
to Q by 

i i i , / Q \ , x 
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FIG. 1: Plots of the universal function F(Q). The line on the left is the small Q approximation in 
(J3.7|) . The line on the right is the large Q approximation in (J3.8|) and (|3.9|) . The square symbols 
are the numerical data of Ref . |40| transformed by (|A9[) and (|A10|) . The plus marks the position of 
the Kosterlitz Thouless transition. 

Note that as Q — > oo, Qd = 3tt/ \n@. For large Q, the expansion of F is 

F(G -> oo) = i- In - | + 0(0* ). (3.9) 

All subsequent terms in the present expansion involve only integer powers of Qjj, with no 
additional logarithms. As discussed in Appendix |XJ the numerical results of Ref. |3^, 
yield the values of F for all values of Q. In particular, at the critical point Q = Q c we have 
from (|A"8|) 

F(£ c ) = 0.502 ± 0.003 (3.10) 

The theory of the Kosterlitz-Thouless transition implies that F(Q) will have a weak essential 
singularity at Q — Q c , similar to that in the specific heat. A plot of the values of F(Q) 
appears in Fig^ It is interesting to note that either the small Q or the small Qd expansions 
is accurate for the entire range of Q values. 

The discussion so far presents our most complete results for the properties of Tql and 
with essentially no approximations. There is, however, still some residual cutoff dependence. 
This can be removed by subtracting corresponding results at two different values of the 
bare coupling R/T (or m*a(T)/T), while U fixed. Depending upon the physical situation, 
changing R/T may also involve some changes in the values of U. However, such changes are 
expected to be small and we neglect any T dependence in U and A in the remainder of this 
section. This allows to obtain an explicit relation between the dimensionless number Q used 
in the present section, and the number g in (|2.2|) . Dividing ()3.3|) by ksT and subtracting 
the corresponding equation at the critical point, and using the definitions in ()3.2j) and (|3.4j) . 
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we obtain 

6 6 1 (TQ C \ 

As expected, g extends from +00 to —00 as Q extends from to +00. Applying the same 
procedure to (jH.fi j) we obtain the universal function in (j2.3j) 

V(g, T/T c ) = A- In + F{G) ~ F{Q e ) (3.12) 



2tt \TQ, 

The expressions (j3.5j) . (j3.7j) - (j3.12j) constitute the central results of this paper. Using as input 
the values of g and T/T c , we compute Q from (j3.11jl and Qd from (j3.8j) : then using results 
(13. 7 j) and (j3.9j) we can compute F(Q), and finally insert in (j3.12j) to obtain T>(g,T/T c ). In 
particular, the small Q expansion in (j3.7j) yields (j2.4j) . Of course, it is better to numerically 
solve for Q from (j3.11j) . rather than obtaining the solution order- by-order in 1/g as was done 
for (gUp . 

We now present some numerical results for the parameters used by Mukerjee and Huse 
@. They set a(T) = a (T - T C MF ). Inserting this in (Q yields 

,= ?#(l-£V (3-13) 



The parameterization a(T) = ao(T — T C MF ) is chosen to be valid near T c , but can also 
be reasonably extended as T — > (in BCS theory, we expect a divergent a(T — ► 0) ~ 
— ln(l/T), but this divergence is expected to be cutoff near a superfluid-insulator transition). 
By p s (0) = -h 2 a(0)/(m*b) in (l3~T3l . we mean the value of the helicity modulus of Zq L 
extrapolated to T = in this manner (the London penetration depth is related to the 
helicity modulus by Xj 2 = l6ire 2 p s (T)/(h 2 c 2 )). It is worth noting here that pJO) and 
T c are, in general, independent of each other, and the Nelson-Kosterlitz relation [38| only 
constrains p s (T c )/T c = 2/n. 

This framework now predicts all ph ysic al properties with 2 input parameters: the values 
of p s (0) and T c . Mukerjee and Huse [35] also defined a parameter 77 as a measure of the 
strength of fluctuations. This is related to the parameters used here by 77 = 2k B T^ IF / p s (0). 
In our numerical results below, we set p s (0) / (k B T c ) = (2/?7)(T c MF /T c ) = 6.8, following their 
parameters. 

An important subtlety should be noted here. The use of (j3.13j) in (j3.11j) normally yields 
a value for Q which decreases monotonically with increasing T, as seems reasonable, given 
our understanding of physical properties of the continuum theory. However, because the 
value of g in (j3.13j) saturates as T — > 00 and because of the presence of the ln(T/T c ) term 
on the right-hand-side of (j3.11j) . for very T the value of Q eventually starts increasing with 
increasing T. This is clearly unphysical, and is an indication that the present continuum 
theory breaks down at large enough T. For the value of rj being used here, this unphysical 
non-monotonicity arises only at T/T c > 20, and we will therefore restrict our attention to 
values of T below this. 

Solving (j3.13j) and (j3.11j) for Q as a function of T/T c , we use the results of this section 
and Appendix |X] to obtain the plot of Fig El for the quantity appearing in (j2.fij) . Note, again 
that either the small Q or the large Q expansion is reasonable accurate. 
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FIG. 2: Plots of the universal function (T/T c )D(g,T/T c ) as a function of T/T c for p s (0) j (&bT c ) = 
6.8. ^Frorn Q we see that (|^| 2 > T = (T/T c )(\^\ 2 ) Tc + (m*k B T c /h 2 )(T/T c )V(g, T/T c ); so (|lf) T 
is determined from the above plot up to an additive, non-singular, linear dependence on T deter- 
mined by (|\I'| 2 )t c - This linear T dependence can compensate for the the linear T dependence in 
the plot above so that (|'I / | 2 )t saturates at high T. Also, as noted in the text, the present theory 
breaks down at large enough T, and its main utility is in capturing the singular increase in (|^| 2 )y 
as T crosses T c . The solid line is the small Q approximation obtained by solving (|3.5|) . (|3,7|) . (|3.10f) . 
()3.11|) . (|3.12|) . and ()3.13|) . The dashed line is the large Q approximation obtained by solving (|3.5I) . 
(EH), (EH), lETTUlh (ETJJ), and KT1 . The square symbols are the numerical data of Ref~E3 

processed via (jlj.llj) . (|3.12f) . (|3.13|) . ()A9j) . and ()A10j) . The plus marks the position of the Kosterlitz 
Thouless transition. 



IV. CONCLUSIONS 

We conclude this paper by discussing some of the experimental and broader implications 
of our work. 

Our primary result ()2.6|) for the coherence length of the charge order can be tested against 
neutron scattering and scanning tunnelling spectroscopy (STS) experiments. However, the 
strong random field disorder may make £^ inaccessible to a neutron probe which averages 
over the entire sample. In contrast, STS provides a local probe, and so may be more sensitive 
to the effects discussed here. 

Consider an STS experiment with a field of view of area A, such as those performed in 
Refs. F l lip. E3. li^. Q uasiparticle interference contributions, such as those computed 



in Ref. [18|, |44|, |45|, |46|, |47|, appear at low temperatures, but we can expect that these will 
significantly broaden at temperartures above T c . We therefore focus here only on the contri- 
bution of the <p fluctuations, which also lead to modulations in the local density of the states 
measured in STS, as shown in Ref. 3 47 ■ We know that the STS measurements are in the 
linear response regime. So, when we perform the Fourier transform of the local density of 
states at the ordering wavevector K c , we find that the signal is proportional to the uniform 
part of the charge order parameter, <p = 0(q = 0). Let us estimate <fi. We have for the free 
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energy 



(FU = , A ^-(hY ,2 U 1,2 + C^ 2 A (4.1) 

Here we used that the result that the random field energy scales as the square root of the 
area. We can now minimize ()4.1|) with respect to <fi 

~ A 1 / 2 ( } 

Taking £|(T) from Eq. ()2.6|) obtain get the temperature dependence of the STS signal at 
the wavevector K c . 

For the case of competition between the superconducting (ty) and charge (0) orders, the 
coupling A in (|2.6|) will be positive. In this situation we have a seemingly counterintuitive 
effect: as T is increased through T c , the amplitude of the charge order is enhanced. The 
physical origin of this is not difficult to understand: the increase in phase coherence as T is 
lowered is associated with an enhanced coherent motion of the Cooper pairs, and this leads 
to a decrease in the amplitude of the spatial modulations [48[. 

An alternative statement of the same physics can be made in terms of the vortices. As 
we argued in Ref. [3, vortices nucleate static charge order, and this was proposed as an 
explanation of the experiments of Ref. (other approaches 0, 0, H^, |2(| 27] have 



proposed static spin order in the vortices — in our theory, static spin order is not nucleated 
by vortices and appears only in phases with global magnetic order (l6|). Increasing T above 
T c causes a proliferation of vortices, and hence an enhancement of charge order. 

While our discussion in this paper has been entirely at the level of the Landau theory 
of multiple order parameters, it is important to keep in mind that such a theory is an 
effective model, and does not preclude other interpretations which focus directly on the 
electronic quasiparticles. In particular we can view the competition between charge order 
and superconductivity as the competition for the ordering of low energy quasiparticles near 
the Fermi surface. So as the superconducting pairing of these quasiparticles is reduced above 
T c , they are more susceptible to charge ordering. 

An interesting direction for future work is to combine the continuum theory of the 
Ginzburg-Landau functional presented here with the theory of time-dependent supercon- 
ducting fluctuations presented in Refs. |3^, 13 Ull this has the prospect of placing more 
precise quantitative constraints on the analysis of the Nernst effect experiments. Moreover, 
the accuracy of either the small Q or large Q expansions suggests that useful analytic results 
may be possible. Results for the fluctuation conductivity in such an approach, inclu ding 
corrections to the Aslamazov-Larkin fluctuation conductivity, have appeared recently [49]. 
A similar dynamic approach can also be applied to computing the linewidths of the elec- 
tronic quasiparticles in the pseudogap regime: the strong amplitude fluctuations in \l/ should 
lead to significant broadening in the electronic spectral functions measured in photoemission 
experiments. 
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APPENDIX A: CORRESPONDENCE WITH THE DILUTE BOSE GAS 



This appendix discusses the connection between the analysis of the dilute Bose gas in 



Refs. |39|, |4(J and the results of Ref. |30| and the present paper. Let us make it clear at the 
outset that we are not advocating a dilute Bose gas description of the underdoped cuprates; 
rather, the finite temperature properties of the dilute Bose gas are characterized by some 
universal numbers which appear also in the models of interest in the present paper. 
The dilute Bose gas is defined by the partition function 



Z B = J T>ip(r, T)e- SB/h 

rh/k B T 



or 2m 2 



(Al) 



We follow the notation of Refs. |3jj, |4Q| throughout this appendix. The only exception is that 
the boson interaction U has been replaced by U B , to prevent confusion with the coupling U 
in ipnjl . 

Integrating out the non-zero Matsubara frequency modes in the Bose gas, the action for 
the zero frequency modes takes the form (|3.1|) with the coupling constants 



h 2 R 

2m 
U 



fi-2U B 
12m 2 U B 



4tt 2 



k B T 



e (h 2 ky(2fn)~^)/(k B T) _ i 



h 2 k 2 /(2 



m) 



(A2) 



The integral above is divergent in the ultraviolet, but if we use ()3.3j) to obtain the value of 
the renormalized coupling R we obtain a convergent integral: 



R + 



2mfi 



4mU B f d 2 k 



ft 



4vr 2 



e {h 2 k^/{2m)- tJ )/{k B T) 



1 h 2 k 2 /(2m)-n 



+ 



2mk B T/h 2 
k 2 + R 



2m 2 k B TU B 



nh 4 



2m 2 k B TU B 

^h 4 



In 



In 



2mfi 



n 2 R^/{k B T) _ V/ 

(2mk B T\ 



\ h 2 R 



(A3) 



In the last expression we have expanded to leading order in U B , as required from consistency 
with previous approximations. Now using the definition of the dimensionless coupling Q in 
()3.4|) . we obtain the value of the chemical potential at the Kosterlitz Thouless transition 



mk B TU B 
nh 2 



In 



mU, 



B 



(A4) 



with universal number £ M computed in Ref. |39| related to the universal number Q c computed 
earlier in Ref. |3(] by 

t> = (A5) 
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Refs. |39j, |40j obtained £ M = 13.3 ± 0.4, which is in reasonable agreement with the value 
g c « 102 obtained in Ref-H the latter value of Q c yields ^ « 14.1 from (lA"5j) . 

The same method can be used to compute the boson density n. Integrating out the 
non-zero frequency modes and mapping onto the classical theory (J3.1|) we obtain 

d 2 k ( 1 T \ m 

d 2 k ( 1 T 2mk B T/K 2 \ mk B T 



4tt 2 ^ e (^ 2 /(2m)-„)/T _ j n 2 k 2 /(2m) -// + k 2 + R J + h 2 F ^ 
= rnksT f 2mksT\ ^ p{Q) (A6) 

in the last equation we have made the same simplification as that in the last equation in 
i3j). The result (|A6|) yields the expression obtained in Ref. for the critical density 



with the universal number £ given by 

e = ^e 27rF ^\ (A8) 

The simulations of Refs. 3^, 4(] obtained £ = 380 ±3, and inserting this result in (jA8|) allows 
us to compute F{Q C ). 

Finally, subtracting (jA7j) from the last equation in (|A6|) we obtain 

71 — 71 S 

c - = —\n(g/g c ) + F(g)-F(g c ^ 



rnksT/h 2n 

= A(X). (A9) 



The function X(X) was computed in numerically Ref. |40j, and its argument X can be related 
to our coupling g by (|A3|) . yielding 

X = 4 + J--ilnf|V (A10) 



g g c vr v ^ 

These earlier results for X(X) therefore yield the needed function F(Q) from (|A9J) and (|A10J) 



APPENDIX B: WEAK AND STRONG COUPLING EXPANSIONS 

This appendix presents discusses the expansion for the universal function F(Q) appearing 
in (|3.fij) for small and large g 

For small g, a simple Feynman graph expansion of ()3.1j) can be carried out to order 
g 2 , with the diagrams shown in Fig El AH propagators in Fig |H] involve the bare 'mass' 
R. A simple calculation shows that the graphs (a,b,c,e) are all absorbed by the first term 
on the right-hand-side of ()3.6|) . after substitution of ()3.3|) . This is a simple example of our 
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(c) 



(d) 



(e) 



FIG. 3: Feynman graph expansion of (|3.1|) for the correlator (|3.6|) . 



claim that all 'tadpole' diagrams can be neglected after substituting R for R. Only FigOJi 
contributes to F(Q) in ()3.6|) . and yields 



F{Q) 



2Q 2 fd 2 p f d 2 q fd 2 k 



9 
9 



4tt 2 
d 2 p 
4^2 



1 



4tt 2 J An 2 (k 2 + l) 2 (q 2 + l)((p + A;) 2 + l)((p + g) 2 + 1) 
1 . / a/p 2 + 4+p" 



2-rrpy/p 2 + 4 



In 



x 



47rp(p 2 + 4) 



v / p 2 + 4 -py 
P+ 2 



y/p 2 + 4 + p 

VP 2 + 4 1 + 4 - p . 



In 



^Bl) 



We evaluated the last integral numerically and obtained (|3.7I) . 

For large Q, the method described in Appendix C of Ref. |30| was used. To one loop order, 
the result 

y ! k B TU 2tt \R d J y 1 

is easily obtained, where, as in Ref.llQ, RdQd = RQ- In obtaining (jB2jL we have to explicitly 
account for all tadpole graphs, and the relationship in (2.4) of Ref. |30| between R and Rb- 
The large momentum behavior of the expansion about (p a = and <p a ^ saddle points 
should be the same, and this ensures that the ultraviolet divergences cancel. At two loop 
order, 35 Feynman graphs appear; these were evaluated as in Appendix C of Ref. 0, and 
their sum was found to vanish. Consequently, there is no order Qu term in F, and the result 

QUID follows. 
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